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' Due to the dispersion of optical phonons, long range electron-phonon correlations renormalize 

\ downwards the coupling strength in the Holstein model. We evaluate the size of this effect both in 

■ a linear chain and in a square lattice for a time averaged e-ph potential, where the time variable is 
5^ ' introduced according to the Matsubara formalism. Mapping the Ifolstein Ifamiltonian onto the time 
^1^, scale we derive the perturbing source current which appears to be non time retarded. This property 

■ permits to disentangle phonon and electron coordinates in the general path integral for an electron 
coupled to dispersive phonons. While the phonon paths can be integrated out analytically, the 
electron path integrations have to be done numerically. The equilibrium thermodynamic properties 
of the model are thus obtained as a function of the electron hopping value and of the phonon 
spectrum parameters. We derive the e-ph corrections to the phonon free energy and show that its 
temperature derivatives do not depend on the e-ph effective coupling hence, the Holstein phonon 

^ ' heat capacity is strictly harmonic. A significant upturn in the low temperature total heat capacity 

over T ratio is attributed to the electron hopping which largely contributes to the action. 

PACS: 31.15.K, 71.38.+i, 63.10.+a 
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I. INTRODUCTION 



, Theoretical investigation on the effects of nonlocality 
Q ■ in electron-phonon coupling [1,2] has grown considerably 
O [ during the last years in conjunction with the large inter- 
est for organic molecular crystals [3], carbon nanotubes 
[4,5] and conducting polymers exhibiting polaronic prop- 
J> . erties. In the Su-SchriefFer-Heeger (SSH) Hamiltonian 
OO ' [6] the electronic hopping is accompanied by a relative 
displacement between adjacent atomic sites that causes 
a nonlocal e-ph coupling with vertex function depend- 
ing both on the electronic and the phononic wave vector. 
This feature leads to peculiar static and dynamical prop- 
[ erties [7] for the SSH model which also displays dimen- 
• sion dependent mass behavior for the acoustical polarons 
' [8,9]. In a path integral investigation of the equihbrium 
, thermodynamics we have shown [10] that the purely elec- 
I ' tronic hopping dominates the low temperature thermo- 
' O \ dynamics as it is responsible for a characteristic upturn 
in the heat capacity over T ratio pointing to a glassy like 
behavior in the one dimensional chain. Moreover, phonon 
anharmonicities induced by e-ph interactions can be large 
in the SSH model due to the time retarded nature of the 
perturbing source current [11]. 
^ In the Holstein model, the coupling of electrons to dis- 

5^ . persionless optical phonons is essentially local and the 
electronic energy at a lattice site linearly depends on the 
atomic displacement at that site [12]. Thus, in the strong 
coupling regime, the unit comprising electron plus in- 
duced lattice deformation is generally a small object (on 
the lattice scale) undergoing a sizeable mass renormaliza- 
tion. Large adiabaticities may however induce a spread 
in the real space for the Holstein polaron which accord- 
ingly becomes lighter. This effect is more pronounced in 
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higher dimensionahty [13]. On the other hand, the dis- 
persive nature of the phonon spectrum is recognized as 
a fundamental feature of the Holstein model as it allows 
for finite values of the electronic site jump probability 
while a dispersionless model would lead to unphysical di- 
vergences [12,14]. Also ground state properties such as 
electron bandwidth and effective mass are more prop- 
erly evaluated versus dimensionality within a dispersive 
model. In this paper we investigate to which extent the 
realistic assumption of dispersive phonons may induce 
non local correlations which consistently renormalize the 
effective e-ph coupling both in one- and two-dimensions. 
We emphasize that the model and the results hereafter 
presented, being valid for any coupling strength, are inde- 
pendent of the existence of small polarons in the system. 

Unlike the SSH model, the Holstein electron hopping 
does not induce a shift in the atomic displacement. Thus 
the perturbing e-ph current does not depend on the elec- 
tronic paths and, along the time scale, it turns out to be 
intrinsecally local. Applying space-time mapping tech- 
niques [15] we show that this feature allows an elegant 
path integral formulation [16] of the Holstein partition 
function in which electron and dispersive phonon coor- 
dinates appear to be decoupled. Section II outlines the 
Hamiltonian model while the path integral method is de- 
scribed in Section III. The results of our work are re- 
ported on in Section IV and some final remarks are given 
in Section V. 



II. THE DISPERSIVE HOLSTEIN MODEL 

We consider the dimension dependent Holstein Hamil- 
tonian consisting of: i) one electron hopping term, ii) an 
interaction which couples the electronic density (///i) to 
the phonon creation (6|) and annihilation (6i) operators 
at a given site 1, iii) dispersive harmonic optical phonons: 

<hm> 

1 

HP^=Y^^{ci)h\h^ (0.1) 
q 

the first sum is over z nearest neighbors, t is the tight 
binding overlap integral. is the Fourier transform of 

h\ and w(q) is the frequency of the phonon with vector 
momentum q. g is the e-ph coupling in energy units. 

The phonons operators can be written in terms of the 
isotropic displacement field as: 

+ ^' = ^ E V2Ma;(q) ^ exp(iq • (1 - n))w„ (0.2) 
q n 



2 



where M is the atomic mass. Then, the e-ph term in 

(1) transforms as follows: 

^'"""^ = ^ E v/2]lMq)^ exp(iq • (1 - n)) 

q l,n 

(0.3) 

The sum over n spans all n*'' neighbors sites of the 1 lat- 
tice site in any dimensionality. Thus, although the Hol- 
stein Hamiltonian assumes a local e-ph interaction, the 
dispersive nature of the phonon spectrum clearly intro- 
duces e-ph real space correlations which renormalize the 
dimension dependent effective coupling. Fourier trans- 
forming the atomic displacement field and taking the lat- 
tice constant |a| = 1, from cq.(3) we obtain for a linear 
chain and a square lattice respectively: 

= ]^ E ///i E ^/2i^Mi) exp(iq' • 1) 

1 q,q' 
X ■Uq'S'd(q' - q) 

SiD{.q' - = 1 + 2 ^ cos(n(«' - q)) 

n=l 

S'2D(q' - q) = 1 + 2 ^ [cosn(g^ - q^) + cosn{q'y - Qy) 

n=l 

+ 2 E [cos(m(g; - qa,) + n{q'y - qy)j 

m,n=l 

+ cos(m(g; - g^) - niq'y - qy)j 

(0.4) 

While, in principle, the sum over n should cover all the 
N sites in the lattice we introduce the cutoff n* which 
will allow us to monitor the behavior of the coupling as a 
function of the range of the e-ph correlations. Note that 
in ID, the integer n numbers the neighbors shells up to 
n* while in 2D, the n* = 1 term includes the second 
neighbors shell, the sum up to n* = 2 includes the fifth 
neighbors shells, n* = 3 covers the nineth shell and so 
on. Switching off the interatomic forces, ti)(q) — loq, one 
would recover from (4) a local e-ph coupling model with 
Sd = l- As no approximation has been done at this stage 
eqs.(4) are general. 

Taking into account first neighbors intermolecular 
forces in the Holstein Molecular Crystal model, the opti- 
cal phonon spectra are given in ID and 2D respectively 
by [14]: 




wii)(q) = ^o/2 + + \l ^tl^ + w^w?i?(q) +^2(2 + /i(q)) 
5(q) = cosqx + cosqy 
h{q) = 2cos{qx - qy) 

(0.5) 
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where wq and ui are the intra- and inter-molecular 
energies respectively. 

Treating in (4) the phonon coordinates as classical vari- 
ables interacting with quantum mechanical fermion oper- 
ators we assume a semiclassical version of the dispersive 
Holstcin Hamiltonian. This is the main approximation in 
the model which will permit us to derive a time depen- 
dent source current for the general electron path integral. 
Averaging (4) on the electronic ground state we also de- 
fine the e-ph energy per lattice site 



< ^r"' >q= tIt^ V^M^^pq,t.q,5,(q' - q) 



which will be used in the next Section. 

While we are thus neglecting the quantum nature of 
the lattice vibrations, the latter may in itself lead to 
retardation effects in the ground state structure of the 
composite quasiparticle made of one electron plus phonon 
excitations in the adiabatic regime. Moreover, at finite 
temperatures, the quantum lattice fluctuations are ex- 
pected to influence the thermodynamics of the system 
mainly for intermediate values of the e-ph coupling [17]. 

III. THE PATH INTEGRAL METHOD 

Let's apply to the Holstein Hamiltonian space-time 
mapping techniques [15,18] which allow us to write the 
general path integral for one electron particle in a bath 
of dispersive phonons. Thus we introduce x(t) and y(T') 
as the electron coordinates at the 1 and m lattice sites 
respectively, and in (1) transforms into 



H%t,t')= -t{f\^{r))f{yir')) + f\y{T'))f{^ir))) 

(0.7) 



subara Green's functions formalism with /? being the in- 
verse temperature hence the electron hops are not con- 
strained to first neighbors sites. Accordingly, eq.(7) is 
more general than H'^ in (1). After setting r' = 0, 
y(0) = 0, we take the thermal averages for the electron 
operators over the ground state of the Hamiltonian thus 
obtaining in d dimensions the average energy per lattice 
site due to electron hopping: 



h%T) ^ ^ ^^>> = -t(G[-x(r),-r] +G[x(r),T]) 




q 



q' 




(0.6) 





(0.8) 
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Vn arc the fcrmionic Matsubara frequencies and ek = 
— 2t ^-^^ y ^ cosiki) is the electron dispersion relation. 

The spatial e-ph correlations contained in (3) are 
mapped onto the time axis introducing the t dependence 
in the displacement field: Uq Uq(T). Assuming peri- 
odic atomic particle paths: Uq(r + /3) = u^{t) we can 
expand Uq(T) in Np Fourier components: 

Nf 



Xqij) = Uo + 2 ({^Un)q cos{LOnT) - (9u„)q sin(w„r) j 

(0.9) 



with u)n = 2mr/(3. 

Then, on the base of eqs.(6) and (9), we identify the 
perturbing source current of the Holstein model with the 
r dependent averaged e-ph Hamiltonian term: 



iir) = E^qM 



iq(r) 



Ar3/2 



v/2ilMq)Euq'(T)pq'5d(q'-q) (0.10) 



q' 



With these premises, we are now in the position to 
write the general path integral [19] for an Holstein elec- 
tron in a bath of dispersive phonons. Assuming a mixed 
representation, the electron paths are taken in real space 
while the phonon paths are in momentum space. The 
electron path integral reads: 



< 



x(/3)|x(0) >= n<x(,3)|x(0) 



>c 



Duci{T)exp - J dT—{u, 



\r) 



+ a;2(q)«^(^)) 



X j D^{T)exp - j^^ dri^^ic^T) + h^r) - jq(T) 

(0.11) 

where the kinetic term (m is the electron mass) is nor- 
malized by the functional measure of integration over the 
electron paths. 

As a direct consequence of the time-local nature of the 
e-ph interactions, the Holstein source current does not 
depend on the electron path coordinates. Then jq(T) in 
(11) can be easily integrated over r using eqs.(9),(10) and 
noticing that 



I 



(0.12) 



we get: 
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Jo 



9d 



(q) = ^V2Ma;(q)^pq.54q' - q) 
q' 



(0.13) 



gd{<l) is thus a time averaged e-ph potential. 

The total partition function can be derived from (11) 
by imposing the closure condition both on the phonons 
(eq.(9)) and on the electron paths, x(/3) = x(0). Using 
(13), we obtain: 



M 



<j) Duq{T)exp 


Puogd{ci) - [ 




Jo 







X j) Dx{T)exp 









m 



x2(r) + /i^(T) 



(0.14) 



Eq.(14) shows that the averaged e-ph coupling is 
weighed only by the r-independent component Uo of the 
displacement field. This property will prove to be essen- 
tial in the next calculations. The integration over the 
phonon paths can be done analytically choosing a mea- 
sure of integration which normalizes the kinetic term in 
the phonon field action: 

V2 



(2Am)(2^^+i) 

Nf 



X JJ (27rn) / di^UnU / d{QUn)c 

n=l 

£'Wq(r)exp|^-— J dTUc^{T) 



(0.15) 



being Am = \J'kTi^(3/M. In the following calculations 

we set M ^ lO^'m. 
Using the result [20]: 



/oo 
dUo exp(6uo — cu^) 
-oo 

we derive: 



exp(6V4c) (0.16) 







j> Dx{T)exp 









^(q) 



exp 



(5d(q)A 



M 



27rw(q)2 



n (2„^)2 + (^^(q))2 



(0.17) 



Eq.(17) represents the main analytical result of our 
model. The exponential function in P(q) embodies the 
effect of the non local correlations due to the disper- 
sive nature of the phonon spectrum. Phonon and elec- 
tron contributions to the partition function are decoupled 
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although the effective potential gd{(l) carries a depen- 
dence on the electron density profile in momentum space 
through the function p(q). 

IV. ELECTRON-PHONON CORRELATIONS 
AND THERMODYNAMICAL RESULTS 

First we analyse the behavior of the time (tempera- 
ture) averaged e-ph potential (eq.(13)) in the case of a 
linear chain and of a square lattice. 

Note that Mq(T) in (9) is a real quantity consistently 
with the closure condition on the lattice displacement 
path. Accordingly, from (6) (and (10)), we take a real 
profile also for the electron distribution in momentum 
space. Setting the total electron density on a single site 
and defining po as the on site electron density, from (6), 
we write the electron profile Pq = po cos(g) for the ID 
system and Pq = Po cos{qx) cos{qy) for the 2D system. 
Since the momentum integration runs over qi G [0, 7r/2], 
Po represents in both cases the total electron density. 
This choice is convenient in order to normalize the av- 
eraged e-ph couplings over the same parameter both in 
ID and 2D. Physically it corresponds to pin one electron 
on a lattice site and to measure the effects on the poten- 
tial due to the e-ph correlations with variable range. 

Let's assume low energy phonon spectra parameters, 
Wo = 20meV and loi — lOmeV. Setting g = 3, wo take 
a strong bare Holstein coupling g although the general 
trend of our results holds for any value of g. In Fig. 1(a), 
giD{q)/po is plotted for three choices of the cutoff n* 
(see eqs.(4)) to emphasize the dependence of the poten- 
tial on the e-ph interaction range. The constant value 
of the potential obtained for cji = is also reported on. 
While in the case of short correlations [n* = 4) there is a 
range of wave vectors in which the effective coupling be- 
comes larger than the dispcrsionless coupling, long range 
e-ph correlations substantially reduce the effective poten- 
tial with respect to the dispcrsionless case and set in an 
oscillating behavior which makes the renormalization q- 
dependent. The potential tends to converge for the value 
n* = 24 which corresponds to 48 lattice sites along the 
chain. 

The projections of the two dimensional e-ph potential 
along the y component of the wave vector is reported on 
in fig. 1(b) for three values of the cutoff on the correlation 
range. For n* = 1 the correlation range is extended to 
the second neighbor shell thus including 8 lattice sites. 
For n* = 2 and n* = 3 we normalize over 24 and 48 
lattice sites respectively. Then, the three values of n* in 
2D span as many lattice sites as the three values of n* in 
ID respectively. This choice permits to normalize consis- 
tently the potential for the linear chain and the square 
lattice. There is a strong renormalization for the 2D ef- 
fective potential with respect to the dispersionless case 
for any value of the wave vector. By extending the corre- 
lation range this tendency becomes more pronounced for 
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Qy close to the center and to the edge of the reduced Bril- 
louin zone. An analogous behavior is found by projecting 
the 2D potential along the axis. The 2D potential sta- 
bilizes by including the 9th neighbor shell (n* = 3) in the 
correlation range. 

Then, an increased range for the e-ph correlations leads 
to an effect which is qualitatively analogous to that one 
would get by hardening the plionon spectrum: a reduc- 
tion of the effective coupling and the consequent lighten- 
ing of the quasiparticle effective mass [21]. 

Let's see now to which extent the momentum depen- 
dent potential affects the thermodynamics of the system. 
In the total partition function (cq.(14)) the electronic 
and phononic paths have been integrated out separately 
as the source current does not depend on the electronic 
coordinates. The q-phonon contribution to the partition 
function is thus given by P{q) in (17) hence, the free 
energy due to the q-mode is: 



where the first two terms represent the harmonic free 
energy f'*(q). Since F''(q) = /S"! ln[2sinh(/3w(q)/2)] 
we find a constraint which allows us to determine the 
cutoff Np = NF{co{q),T). For the ID system, at 
T = IK and for the phonon parameters given above, 
we get Np ^ clO^ with c varying in the range [4.2 — 4.5] 
according to the g-mode. 

The e-ph interactions renormalize downwards the har- 
monic values for any q as Fp{q) — F'^{q) cx — g"^ (q) / uj"^ (q) 
but the temperature derivatives of the free energy do not 
involve the coupling term. Then, on general grounds, 
phonon entropy and heat capacity are not affected by 
the strength of the e-ph interaction and signatures of e- 
ph anharmonicity should not be expected in the equi- 
librium properties of the phonon subsystem. This is a 
direct macroscopic effect of the local nature of the inter- 
actions in the Holstein model. On the other hand lattice 
non linearities may appear in the ground state proper- 
ties of the Holstein Hamiltonian and significantly modify 
the dynamical properties of the polaronic quasiparticles 



In Fig. 2, we consider the linear chain with the cut- 
off n* = 24 of Fig. 1(a) and set po = 1- The phonon 
free energies are plotted for two selected wave vectors 
together with the corresponding harmonic values. It is 
seen that the renormalization of the q-dependent free en- 
ergies due to the e-ph coupling is very weak. However, 
being 5^(q)/ti)^(q) cx g^/a)(q), the phonon free energy 
reduction may be appreciable in very low energy phonon 
spectra provided that the bare coupling g is sufficiently 
strong. 

Next we turn to the evaluation of the electronic con- 
tribution to the total partition function (eq.(17)). The 




(grf(q)^)^ 

2Mw(q)2 



(0.18) 



[22,23]. 



8 



results hereafter presented refer to the case of ID elec- 
tronic paths. As the hopping energy density is related 
to the electron propagator (eq.(8)) a preliminar wave 
vector integration prepares as a function of the path 
coordinates. Then, the § Dx{t) integration has to be 
done numerically. Assuming the periodicity condition 
x{t) = x(t + (3), the particle paths can be expanded in 
Np Fourier components 

x{t) Xo + 2 {^Xn cos(^„r) - 9a;„ sin(^„r)^ 

n=l 

e„ = 2^n//3 (0.19) 
so that the functional measure of integration 

{2Xmr^'-+'^> J-A J-2A J-2A 

Xm = \J'Kn^l3/m (0.20) 
normalizes the free electron term in (17): 



^ Dx{T)exp — J dT^x.'^{T 



1 (0.21) 



Eq.(21) provides a criterion to set the cutoff A which 
turns out to be cx Xm thus implying that the maximum 
amplitude of the particle path coefficients becomes large 
in the low temperature range. This is consistent with 
the physical expectation that towards high temperatures 
only a reduced number of paths contributes to the ac- 
tion which thus tends to the classical limit. Accordingly, 
the integration over each Fourier component in (20) is 
carried out over a number of points A^a cx K/Vt and 
the electron hopping in the path integral (18) is weighed 
over a total number of electron paths (A''a + 1)2^p+i. We 
find that K ^ 20 and Np = 2 suffice to ensure numerical 
convergence. While the free electron term only scarcely 
contributes to the electronic action, the hopping depen- 
dent action is dominant also at low temperatures. Infact, 
according to our measure of integration, the ensemble of 
relevant particle paths (over which the hopping energy 
density is evaluated) is T dependent. However, given a 
single set of path parameters one can monitor the be- 
havior versus T. It turns out that the hopping decreases 
by decreasing T but its value remains appreciable also at 
low temperatures. Since the dr integration range is larger 
at lower temperatures, the overall hopping contribution 
to the total action is relevant also at low T. 

The equilibrium thermodynamics of the Holstein 
model can be derived from (17). We average the phonon 
contribution over a set of 50 modes in the ID Brillouin 
zone and take the bare e-ph coupling g = 3. Phonon 
and electron free energies are plotted separately versus 
temperature in Fig. 3(a). Note that in the low T range 
the electron free energy has a positive derivative both for 
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narrow and wide band values although the range shrinks 
by increasing the hopping integral value. For t — 0.2eV, 
the electron free energy gets the maximum at T ^ 60K 
while for t = AeV, the maximum occurs at T ~ 13K. 
This feature is mirrored in the total heat capacity which 
is reported on in Fig. 3(b) for four t values together with 
the phonon heat capacity normalized over the number 
of modes. As noted above the latter overlaps the purely 
harmonic heat capacity no matter how large the e-ph cou- 
pling may be. Thus the Holstein phonon heat capacity 
does not contain signatures of e-ph anharmonicity. The 
total heat capacity over T ratio shows an upturn at low 
temperatures as shown in Fig. 3(c) due to the electron 
hopping mechanism which has already been encountered 
in the study of the SSH model [11]. However, while in 
the latter the shape of the low T upturn could be tuned 
by the strength of the e-ph coupling in the source ac- 
tion, in the Holstein model the low T anomaly has purely 
electronic origins as the e-ph effects are contained in the 
phonon partition function. This also explains why the 
absolute values of the Holstein heat capacity are much 
lower than the SSH heat capacity throughout the whole 
temperature range. 

V. CONCLUSIONS 

In the Holstein model, the e-ph coupling is local and 

the vertex function does not depend on the wave vectors. 
Nonetheless non local e-ph correlations arise in the sys- 
tem once the dispersive nature of the phonon spectrum is 
taken into account. I have studied a semiclassical version 
of the one electron Holstein Hamiltonian in which quan- 
tum mechanical fermion operators interact with classical 
lattice displacements. Applying the path integral method 
the real space interactions can be mapped onto the time 
scale and, in the Matsubara formalism, we have derived 
an analytical expression for the time averaged wave vec- 
tor dependent e-ph coupling. Assuming a force constant 
intermolecular model we have evaluated the renormaliza- 
tion of the wave vector dependent e-ph coupling, both in 
one- and two-dimensions, as a function of the range of 
the e-ph correlations. In ID, for short range correlations 
(n* = 4), I find an enhancement of the effective coupling 
with respect to the dispcrsionless coupling in a window 
of phonon wave vectors. By increasing n* , the effective 
coupling renormalizes downwards for any q although the 
reduction is larger towards the edge of the Brillouin zone. 
In 2D, the renormalization is even more pronounced and 
takes place also for small values of the cutoff on the cor- 
relation range. 

Thus, the inclusion of the e-ph correlations in the ef- 
fective coupling has an effect which is similar to that 
one would obtain by enhancing the characteristic phonon 
frequency and it ultimately leads to reduce the effective 
mass of the quasiparticle. The Holstein polaron mass 
in a dispersive phonon model is thus lighter than in a 
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dispersionless model. 

Mapping the e-ph interaction onto the time scale we 
also find that the perturbing source current of the model 
depends on the time only through the atomic displace- 
ment field, hence it is not retarded. Then, the temper- 
ature (time) averaged current is proportional to the av- 
eraged atomic displacement and it does not include the 
electronic coordinates. This property allows us to decou- 
ple phonon and electron degrees of freedom in the general 
path integral and in the total partition function. Such 
a disentanglement does not occur in the Su-Schrieffer- 
Heeger model [10] as the source current does depend on 
the electron path coordinates. As a physical consequence 
the renormalized Holstein e-ph coupling is frozen in the 
phonon partition function while the electron action is 
dominated by the hopping energy density. After inte- 
grating out analytically the phonon degrees of freedom 
we have obtained the q-dependent phonon free energy 
and estimated the e-ph corrections (with respect to the 
harmonic values) in the thermodynamical properties. It 
turns out that the Holstein heat capacity does not show 
any signature of phonon anharmonicity induced by the 
e-ph interactions. This marks a striking difference with 
respect to the thermodynamical behavior of the previ- 
ously investigated Su-Schrieffer-Heeger model. On the 
other hand, the electron hopping energy strongly con- 
tributes to the total action of the Holstein model and, 
at low temperatures, it is responsible for a peculiar up- 
turn in the heat capacity over temperature {Cy/T) ratio. 
While qualitatively the latter feature had also been en- 
visaged in the SSH heat capacity, in the Holstein model 
the shape of the broad upturn has purely electronic ori- 
gins and it is not affected by the e-ph source action. This 
explains why Cv /T has a negative temperature deriva- 
tive throughout the whole T range for any value of the 
hopping integral. 
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FIG. 1. (a) Temperature averaged e-ph coupling (in units 
meVA~^ versus wave vector for a linear chain, n* represents 
the cutoff on the e-ph correlations, po is the electron density. 
LJo — 20meV and u>i = lOmeV. g — 3. The dispcrsionloss 
e-ph coupling is obtained for ujo = 0. (b) Temperature aver- 
aged e-ph coupling in two dimensions versus the y component 
of the momentum. The cases n* — 1, n* — 2 and n* = 3 
imply that the correlation range includes the second, the fifth 
and the nineth neighbors shell, respectively. The input par 
rameters are as in (a). 
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FIG. 2. Anharmonic Fp and harmonic F*^ phonon free 
energies for two values of the one dimensional wave vector 
g = 0, 7r/2. The phonon spectrum parameters are as in 
Fig.l(a). 
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FIG. 3. (a) Phonon free energy (symbols) and electron free 
energies for four values of the hopping integral t. The phonon 
energies and the e-ph coupling are as in Figs.l. (b) Total 
(electron plus phonon) heat capacities for four values of t. 
The phonon heat capacity is plotted separately, (c) Total 
heat capacities over temperature ratio. 
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